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We study the effective diflFusion constant of a Brownian particle linearly coupled to a thermally 

^— ( ■ fluctuating scalar field. We use a path integral method to compute the effective diffusion coefficient 

^— ( ' perturbatively to lowest order in the coupling constant. This method can be applied to cases where 

' the field is affected by the particle (an active tracer), and cases where the tracer is passive. Our 

04 . results are applicable to a wide range of physical problems, from a protein diffusing in a membrane 

t I ' to the dispersion of a passive tracer in a random potential, fn the case of passive diffusion in 

^ l' a scalar field, we show that the coupling to the field can, in some cases, speed up the diffusion 

, corresponding to a form of stochastic resonance. Our results on passive diffusion are also confirmed 

■ via a perturbative calculation of the probability density function of the particle in a Fokker-Planck 
CO ' formulation of the problem. Numerical simulations on simplified systems corroborate our results. 

. PACS numbers: 

o . 

CD ■ I. INTRODUCTION 

I K ' Diffusion in a quenched random medium is a problem which has been extensively investigated [l|, 0| • An important 
, physical problem is to understand how the long-time transport properties of a Brownian particle, notably its effective 

■ diffusion constant, are modified with respect to those of a homogeneous medium. Two types of problem with quenched 
^ ' disorder have been extensively studied. The first is diffusion in a medium where the local diffusion constant depends 

on spatial position and is taken to have some statistical distribution and the second is for a particle which is advected 
Q ' by a quenched random velocity field. The case where the random velocity field is derived from a potential is of 
importance as a toy model for spin glasses and glasses as it can exhibit a static spin glass like transition |^ Q and 
a dynamical structural glass-like transition [SflZl]- Variants of the toy model of diffusion in a random potential can 
exhibit a transition in their transport properties notably diffusion which is normal in the high temperature phase and 
becomes subdiffusive in the low temperature phase P, I8l-[l0| . The onset of the subdiffusive regime is signaled by the 
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vanishing of the late time diffusion constant. It should also be mentioned that above class of problems can also be 
related the problem of evaluating the macroscopic/effective electrical properties, such as conductivities or dielectric 
constants, of random conductors or dielectrics 
. The case of diffusion in dynamically evolving random media has perhaps received less attention. The most widely 
' studied problem of a time dependent nature is for a particle diffusing in a turbulent flow [ill . [T2j , and the particle 
here is a passive tracer and has no effect on the flow field. Another example is for a protein diffusing in a membrane, 
04 where the protein is subject a force generated by membrane curvature or composition, and these quantities themselves 
fluctuate. However in this case the membrane is in general affected by the protein [isj . Naively one might expect 
' that protein diffusion is speeded up by coupling to height or composition ffuctuations in a membrane, however this 
is not the case, the feedback of the protein on the membrane configuration actually slows down the protein diffusion 
with respect to a non-fluctuating homogeneous one. A general question that arises in these sorts of problems is: when 
does the fluctuating field speed up the diffusion of a tracer and when does it slow it down ? An important question, 
, , which has received much recent attention, is how the diffusion constant of a protein depends on its size. The classic 
hydrodynamic computation of Saffmann and Delbriick pd| treats the protein as a solid cylinder in an incompressible 
I layer of fluid (the lipid bilayer) sandwiched between another fluid (the water). In this formalism the protein diffusivity 
" ■ ■ shows a weak logarithmic dependence on the cylinder radius. However the validity of this result has been called into 
question experimentally where a stronger dependence is reported [T5| . A number of theoretical studies have suggested 
that Oa protein's diffusion could be modified by coupling to local membrane properties, such as composition and 
curvature, that are not taken into account in a purely hydrodynamic model for a membrane [isl [l6l - l22| . One should 
bare in mind that protein coupling to local geometry and composition has also been postulated as a mechanism for 
inter- protein interactions in biological membranes [Ij, . 

In this paper we will consider a particle, whose position is denoted by x{t), diffusing in a scalar field hip, where h 
is a parameter controlling the coupling of the field to the particle. Thus the particle drift is given by it = — k^^Vi/', 
where denotes the bare diffusion constant of the particle. The dynamics of both the particle and the field are 
overdamped and stochastic. Our aim is to evaluate how the effective diffusion coefficient for the particle is modified 
by its coupling to the field. A path integral approach allows us to perform a perturbative computation if the coupling 
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constant between the particle and the field is small for a wide range of physically different problems. Within this 
formalism we recover some previous results on the diffusion of active tracers where the whole system (particle plus 
field) obeys detailed balance [l^l • However the method also allows us to analyze passive diffusion and a continuum of 
intermediate models with varying degrees of feedback of the particle on the field and for thermal fluctuations on the 
particle and field that are not necessarily at the same temperature. This intermediate range of models could apply to 
the study of tracers in non-thermally driven fields, and have applications to active colloidal systems [l^, swimmers 
[2^ or systems where the tracer is heated with an external heat source such as a laser The basic method may 
also be useful to study systems where the field itself is non-thermally driven, for instance lipid-bilayers where the field 
is driven an electric field [28j . 

The formalism also allows us to see how the relative time scales between the tracer and fluctuating field affect the 
effective diffusion constant of the tracer. 

Our results for the effective diffusion constant show very different effects for active and passive diffusion, and 
diffusion in a quenched potential. In pure active diffusion, where the stochastic equations of motion obey detailed 
balance, the diffusion is always slowed down. For passive diffusion, if the field evolution is slow, the particle slows 
down, as if was evolving in a quenched potential [2!|. As the field evolution is speeded up the diffusion constant 
increase and then reaches a maximum, in a manner reminiscent of stochastic resonance (29ri32] . As the field evolution 
rate increases still further, the diffusion constant diminishes and eventually reaches the bare value it would have in 
the absence of coupling to a field. The results derived via the path integral formalism for this case are also derived 
via perturbation theory on the corresponding time-dependent Fokker-Planck equation. 

Finally we demonstrate the analytically predicted effects by numerical simulation of some simple models where the 
fluctuating field has a small number of Fourier modes. 



II. THE MODEL 



To start with we deflne the general class of models of interest to us in this paper. First we consider the dynamics 
of a Langevin particle whose position is denoted by x{t) diffusing in a d dimensional space with a linear coupling 
to a fluctuating Gaussian field, taking values of the same d dimensional space, as introduced in [l^. The overall 
energy/Hamiltonian for the system is 

i? = i y (l){x)A(j){x)dx - hK(j){x{t)), (1) 

the first term is the quadratic energy of a free scalar field and the second corresponds to the tracer seeing an effective 
potential hip — —hK(f>. We will take A and K to be self adjoint operators. For two operators A{x,y) and B{x,y) 
we will denote by AB{x,y) their composition as operators: {AB){x,y) = J A{x, z)B{z,y)dz. 

For a system obeying detailed balance (and whose equilibrium state is thus given by the Gibbs-Boltzmann distri- 
bution) the particle evolves according to 

SH 

x{t) = -Kx^ 1- VK^vit) = hKr^VK<j){x{t)) + y^r]{t), (2) 

ox 

where the noise term is Gaussian with zero mean zero and correlation function 

{ri{t)rj{sf)^2T5{t~s)l, (3) 

where T is the temperature of the system. 

In the absence of a coupling between the field and the particle (/i = 0), the particle diffuses normally and the mean 
squared displacement at large times behaves as 

({x{t) - x{0)y) - 2dTK^t = 2dD^t, (4) 

where d is the spatial dimension and Dx ~ Tk^ is the bare diffusion constant. 
We take a general dissipative dynamics for the field [33| : 

(j){x, t) = -f ^/K^^(a;, t) = ~K^RA(l){x) + hK^RK{x - x{t)) + ^/k^^{x, t), (5) 

where i? is a self adjoint dynamical operator and ^ is a Gaussian noise of zero mean which is uncorrelated in time. In 
order for the field to equilibrate to the Gibbs-Boltzmann distribution, the correlation function of this noise must be 

(C(a;, my, s)) = 2TR{x - y)S{t - s). (6) 
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The effective diffusion constant for the particle defined via 

({x{t) - X{0)y) - 2dTKet=2dDet, (7) 

De is the effective or late time diffusion constant. 

Our model applies to many systems: a point magnetic field of magnitude h diffusing in a Gaussian ferromagnet 
can be modelled with A{x,y) — (— + m'^)6{x — y), K{x,y) = d{x — y) and R{x,y) = 6{x — y) for model A 
dynamics or R{x,y) — — V^(5(a; — y) for model B conserved dynamics [s^. Here the field (j) can correspond to the 
fluctuations of a range of order parameters in a lipid membrane, in the Gaussian approximation, where the fluctuations 
are weak 20, 21]. The order parameter in question could be compositional fluctuations above the demixing transition 
in lipid bilayers or binary fluids. Here, a delta function form of the coupling K would correspond to the tracer 
having a preference for one of the lipid phases or preferential wetting for one of the phases of a binary fluid. As well 
as compositional fluctuations, </> could also correspond to fluctuations in local lipid ordering (gel/liquid phases) and 
possibly orientational order in the lipid tails. It could also correspond to local membrane thickness in the case where 
the protein induces a local hydrophobic mismatch and locally alters the thickness of the bilayer. 

To model a lipid membrane where the field represents the height fluctuations and the particle is a protein coupled 
to membrane curvature, we may take the Helfrich Hamiltonian ^ A{x,y) = (kV^ - crV^)(5(a; - y), K{x,y) = 
— V^(5(a; — y) and R is given by its Fourier transform, R{k) = l/4ry|fe|, where 77 is the viscosity of the solvent 
surrounding the membrane (ssl . ISq ] . 

The computations below can be made a little more general. As mentioned in the introduction there are a number of 
physical cases of non-equilibrium systems where one is not restricted to stochastic dynamics obeying detailed balance 
and for instance the coupling between the particle and the field may be taken non-symmetric: instead of ([5]), we will 
consider 

^{x, t) = -K,^RA(t){x) + ChK^RK{x - x{t)) + y/^£,ix, t), (8) 

where we have introduced the parameter C. The case of stochastic dynamics with detailed balance is recovered for 
^ = 1. We note that Eq. ([2]) has been extensively studied in the case where the field (p evolves independently of the 
particle position, which corresponds to ^ = 0. This problem is referred to as the advection diffusion of a passive scalar 
(the particle concentration) in a fluctuating field (j). It had been suggested that this form can be used to approximate 
the diffusion of an active tracer particle in [l6l [l8| , for the case of a protein weakly coupled to membrane curvature. 
In this approximation it was found that the effect of the field fiuctuations could be to increase the diffusivity of the 
tracer particle with respect to that obtained when it is not coupled to the fiuctuating field {h — 0). However the 
numerical simulations of [l^ where the effect of the particle position on the field is taken into account showed that 
the diffusion is reduced with respect to the case h — 0, in agreement with later analytical studies (lol [2^ 

Also as previously mentioned, a further generalization can be made to our model: since, in a non-equilibrium 
system, the particle and the field are not necessarily driven by the same thermal bath, they can experience different 
temperatures. In this case he temperatures appearing in the correlation function of the noises in Eqs. ([3|) and ^ will 
be denoted respectively and T^. 



III. EFFECTIVE DIFFUSION EQUATION AND PATH INTEGRAL FORMALISM 



A. Effective diffusion equation 



Our aim is to study the average value of the mean squared displacement of the particle, we thus integrate the 
dynamical equation of field Eq. assuming without loss of generality that the field = at time t = 0. This gives 



(t>{x,t)^ e-''i-^*-'^"^^[ChK^RK{x-x{s)) + ^£_{x,s)]ds. (9) 

>.' — 00 

Using this result in Eq. ([2]) we obtain the effective diffusion equation for the particle: 

x{t) = hK^VK I e-''*(*-'^)'RA \j;^hK^JiK{x(t) ~ x{s)) + ^i{x(t), s)] ds + y/J^r]{t). (10) 

The right hand term can be split into two parts: a deterministic part depending only on the particle trajectory, 
and a stochastic part which depends on the noise driving the field and on the particle trajectory: 



x{t) = y^r]{t) + f F 

J —00 



{x{t) - x{s), t - s)ds + a{x{t), t), (11) 



4 



where H is a Gaussian noise dependent on the position in space and time, with correlation function 

{Six, t)H(y, sf) = TG{x -y,t~s), (12) 

and we have introduced the functions 

F{x,u) = Ch^K^K^VKe-''^'^"^RK{x), (13) 
G{x,u) = -h^nlVV'^K^e-''*'\''\^^A-\x). (14) 

We introduce the above functions for two reasons: provide more compact notation for the switch to the path 
integral formalism, and show explicitly which part of the following computation is general and does not depend on the 
expressions of F and G. Indeed, for a different choice of functions, dTTI) and (|T2|) define a more general model, which 
may be analyzed using the method presented here. We will need explicit expressions for F and G; in particular their 
X dependence is rather obscure. Fourier transforming allows us to write them as a sum of functions with a completely 
explicit X dependence: 

F{x,u) = Ch^n^K^ I |l^»fce-«^"^('=)^(^)E(fe)if2(fc)e''^- = J ^FUx,u) (15) 



and 



G{x,u) ^h^Ki I ——rkk^ = \__!_^thx^ I _^^Q \ 

^ ^ " ' {2Tr)d A(fc) J {2TT)'i '"^ ' ' ^ ^ 



B. Path integral formulation 

We now turn to the path-integral formalism; the first steps are analogous to those described for general stochastic 
dynamics in [37l - l39l | and [40l | for transport by a time dependent incompressible velocity field. The partition function 
for this system is: 

Z = J (^x{t) - y^Tjit) ~ J F{x{t)-x{s),t-s)ds-S{x{t),t)^ P[ri]Q[S.][dx][dTi][dS], (17) 

where P and Q are the functional Gaussian weight for the noises. Note that we should in principle include a Jacobian 
in the above expression. However if we use the Ito convention the Jacobian term is constant and independent of 
the path. This is because the causality of the equation of motion and the use of the Ito calculus mean that the 
transformation matrix is triangular with diagonal terms which are constant (391] . We then use an functional integral 
description of the S function, introducing the vector field p, 



Z = y exp y p(i) • x{t) - y/K;;:r]{t) - J F{x{t)-x{s),t- s)ds~S{x{t),t) 



dt] P[rj]Q[S][dx][dp][drj][dS] 



(18) 

Then we use the standard result that, for a Gaussian random variable of zero mean u, (exp(aw)) = exp(a^ (m^)/2) to 
perform the integration over the noises to obtain 



Z ^ J cxp{-S[x,p])[dx][dp] (19) 

J F{x{t) ~ x{s),t- s)ds^ dt+Dx J \p{t)\'^ dt+'^ J p'^ {t)G{x{t)-x{s),t-s)p{s)dt ds. 

(20) 

This action is the sum of the action of the pure Brownian motion 

So[x,p]=-i j p{t)-x{t)dt + D^ j \p{t)\^dt (21) 



where the action is given by 
S[x,p\^-i I p{t)-(x{t)- 
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and the action of the interaction 



Sint[x,p] = i J p{t) ■ F{x{t) - x{s),t- s)9{t- s)dtds+'^ J p^{t)G{x{t) - x{s),t - s)p{s)dtds. (22) 

Since G{—x,—u) = G{x,u), we can write this integral only with times satisiying t> s, which will be convenient for 
the ensuing calculations: 

,p]=i j p{t)-F{x{t)-x{s),t-s)e{t-s)dtds + T^ j p'^{t)G{x{t) - x{s),t - s)p{s)0{t - s)dtds. (23) 



C. Computing averages with the free action 

We will need to compute averages with the free action Sq. Since it is quadratic in x and p, we just need the one 
point and two point correlation functions. Moreover, the position of the particle is relevant only with respect to its 
position at, say, t = 0. In order to keep compact notations we define 

xo{t) = x{t) - x{Q). (24) 

The correlation functions required are: (£Co(i))oj {p{i))oi {p{t)p{s)^)o, {xo{t)p{s)'^)o and {xo{t)xo{s)^)o. 
Using the symmetry of ^o, we have immediately 

{xo{t))o = 0, (25) 
(p(0>o = 0. (26) 

To obtain the two point correlation functions, we use the fact that the (functional) integral of a total (functional) 
derivative is zero, for example 

= 1 ^ (p(t)e-^°) [dx][dp] = J p{t)p{sfe-^'^[dx][dp], (27) 

dividing each side by Zq, we obtain (p(t)p{s)^)Q = 0: {p{t)p{s)^) is a constant. This constant must be zero, 
because the action So does not correlate p at different times: the p{t) are all independent. We thus have the first 
correlator 

{pit)pisf)o = 0. (28) 
We use the same technique for the other correlators: 

= 1 (p(t)e-^'') [dx][dp] = J {S{t - s)l+p{t) [ix{sf - 2DMsf])e-'''[dx][dp], (29) 

which gives, after dividing by Zq and integrating over s G [0,t], 

{xo{t)p{sf)^ = 1X10,4^)1, (30) 

where Xa{s) is the characteristic function of the set A (equal to 1 if the argument is in A and zero elsewhere). Finally 
the identity, 

= 1 {xo{t)e-^°) [dx][dp] = J xo{t) [ix{sf - 2DMsf] e-^°[dx][dp], (31) 

leads to 

{xo{t)xoisf)o = 2D^L{[0, t] n [0, s])l, (32) 

where L{I) is the length of the interval / C M, this is of course the standard result for free Brownian motion 

(^Xo{t)xo{s)'^)^ = 2Dx mm{t, s) ii t,s > 0. 
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IV. PERTURBATIVE CALCULATION OF THE EFFECTIVE DIFFUSION CONSTANT 



A. Derivation of the general result 

Here we should go back to our aim of computing the effective diffusion constant Df.. To do this, we have to evaluate 
(a;o(t/)^) at a large time tf. Since we do not know how to compute averages with the action S, we use a perturbation 
expansion in terms of averages over 5*0, which will be denoted (. . . )o: 

(cxp(-S'i„t[a;,p]))o 

Averages appearing here are not easy to compute, but it is easily to compute the first nontrivial term in expansion in 
the interaction action S'lnt- To do this, we just expand the exponential functions: 

^""^'^■^ ^ (l-^i„t[.,p]>o ■ ^''^ 

The interaction action is linear in F and G and so are the averages in ([Ml) ; these functions are the sum over Fourier 
modes of functions Fk and Gfc, so we can carry out the computation with only one Fourier mode, and integrate over 
all modes at the end. We note S'lnt, fc the interaction action associated to the fc-mode, and compute (Sint.fci^c, p])o and 
(a;o(i/)^S'i„t^fc[x,p])o. 

In what follows we work at fixed wave-vector fc, so A(fc), K{k) and R{k.) are pure numbers and to lighten the 
notation we will write them A, K and R. 

Every average we have to compute is made of terms of the form (IljLi Oj-e*'"'^), where Oj are operators linear in 
X and p. We will need the following formula, which is easy to derive from the Wick's theorem, 

(n o.e^'A = e-'-^^(-^)^ E f n ■ (o.-> (n o.) V (35) 

\i=i / JdN \jeJ \]<^j I j 

where N is the set {1, . . . , n} and the sum over J denotes the sum over all subsets of TV. 

We start with (S'i„t.fc[a;,p])o- This average contains two integrated terms. The first is • Fk{x{t) — x{s),t — s))q 
with t > s, which involves 

^p(t)e''=-("(*)-"(^»)^ =2(p(0(a;(t)-a?(s))^)Qfce-'='^-^*-^) =0, (36) 

with the notation fc^ — |fcp and we have used that {p{t)x{t)'^)Q — because we use the Ito convention in our path 
integral. For the same reason, (j){t)^ Gk[x{t) — x{s),t — s)p{s)'^^ = 0. Hence 

{SintAx,p])o = 0. (37) 

Now we turn to {xo{tf)'^Sint,k[x,p])o- We have to compute (a;o(i/)^p(i)e*'"''^^'^*^~"''^^)))Q and 

(a;o(i/)^p(i)p(s)"^e*'"'^^'^*^^"''^*^-')p, with s < t. The only non-zero term in the first average is 

{xoitf)Mty''-^^^'^-^^'^^)o = 2t{p{t)xo{tf)^)o{xo{tf)ix{t) - a;(s))^)ofce-'='^^(*-«) 

= ~AD^L{[0,tf]n[s,t])ke-'''^^('-^)xio,tf]{t) (38) 
= -4i?,(t-max(s,0))fce-^-'^-(*-^)x[o,t,](i)- 

The second contains two non-zero terms and we get 

(a^o(i/)'p(t)p(s)^e*'=-("W-="(^»)^ = [4i?.(t-max(s,0))fcfc^-2x[o.t,](s)l] e-'^-'^-(*-^)x[o.t,] (0- (39) 

Now we just have to integrate the above results, and since we are interested in the long time behavior, we can 
neglect the terms in o{tf): 



ixQ{tf f j p{t) ■ Fk{x{t) - x{s),t- s)e{t - s)dtds^ = 
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and 



^xoitjr J P^{t)G^{x{t) - x{s),t- s)pis)dtds'j^ = 2T,h\lA-'K^^^^^^^^^tf. (41) 



Thus, gathering these two resuhs in p4| gives 



(xoitf) ) ^ 2dtf [D^-hn^kK ) . (42) 



Integrating over the modes, we get for the effective diffusion constant 

^jj _']^ fc2j>.^ ^2 ^xT^£'c.fc' + - A^.r^)^0i?(fc)A(fc) 

" rfi (271)'^'"" ^ ' A(fc)[i?,fc2 + K^i?(fc)A(fc)]2 ■ ^ ^ 

This expression is our main resuh. The large number of parameters makes its interpretation quite difhcuh, so we will 
apply it to some special cases to show the great variety of phenomenon which could be described. 



B. Application to some special cases 

1. Stochastic dynamics with detailed balance 

We first analyze the case where the particle and the field see the same temperature, — = T, and the dynamics 
obeys detailed balance, i.e. C = 1- The effective diffusion constant is thus 

D'^^ = D (l-l^f ^ "-fe'^(fe)' ^ (44) 

This result is exactly what was found in j22l] via a Kubo formula formalism which we emphasize applies only to this 
particular case. By inspecting Eq. (pi)) we see that the correction to the bare diffusion constant is always negative, 
the diffusion is thus slowed down by its coupling to the field. The fact that the diffusion is slowed down for all values 
of h, and not just in the regime of small h can be shown explicitly within the Kubo formalism [2^ . A physical 
explanation for the slowing down of diffusion in this case can be found in studies of the drag on a particle coupled to 
a field. The reaction of the field to the particle is to create a polaron like deformation of the field about the particle, 
however a moving particle has a polaron which is not symmetric with respect to the front and rear of the particle. 
This deformation generates a drag force which tends to pull the particle backwards [lO, [HI . 



2. Passive diffuston 

For passive diffusion, that is C = 0, and still with equal temperatures for the particle and field , we have 

^pa. ^ ^ ^« e Kit? - ^^^(fe)^(fe) ^ (45) 

d J {2nY ^ ' A{k)[Dxk^ + K^R{k)A{k)Yj' ^ ' 

This result shows that, depending on the speeds of evolution Kx and of the particle and the field, the particle may 
speed up or slow down. This is in agreement with the following intuition: with a slow field we are close to diffusion in 
a quenched potential, which slows down the diffusion due to trapping in local minima which are temporally persistent. 
However when the field fiuctuates quickly the field fluctuations kick the particle along thus adding to the effective 
random force the particle experiences. However this picture is not totally valid, there appears an optimal value of 
at which the perturbative enhancement of the tracer's diffusion constant is maximal. If the field fluctuates too quickly 
the effect of field fiuctuations simply average out to zero. This result is quite difficult to understand physically but it 
resembles closely the phenomenon of stochastic resonance [29l - [3l| , where the application of a periodic but deterministic 
potential to Brownian particles can show an optimal frequency at which the particle dispersion is maximized. Here, 
the optimal value of depends, in particular, on the bare diffusivity of the particle, thus two different species 
will react differently to the same external field. We note that this type of phenomena can be used to sort molecules 
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3. Particle not connected to a thermal bath 



Another case which could have physical relevance is where the particle is not connected to a thermal bath, i.e. 
— and energy enters the system only via the fluctuations of the field. Hence in the absence of coupling to the 
field, the particle cannot diffuse, as ~ 0. Our result shows that the fluctuations of the field will induced a non-zero 
diffusion constant for particle: 

The effective diffusion constant is now proportional the temperature seen by the field T^: in some sense, the field 
acts as a thermal bath for the particle. Note the effect of field fiuctuations is to speed up the diffusion from a zero 
diffusion constant to that given by Eq. ()46|) . this is in agreement with the physical intuition that the field fluctuations 
will help the particle to disperse. Interestingly we see that when = the result for is independent of C, and 
the active and passive cases have the same diffusion constant. 



V. NUMERICAL SIMULATIONS 



In this section we test out our theoretical predictions against the numerical simulations of a toy model in one 
dimension. We should bare in mind that the simulation of the diffusion of active tracers is more computationally 
intensive than that for passive tracers. In the latter case we can simulate the diffusion of an ensemble of independent 
tracers on a given dynamical realization of the fluctuating field. In the former case however we must follow the 
diffusion of a single particle in the fluctuating fleld as for a system with more than one particle the coupling to the 
field introduces interactions between the particles [H, 111] . 



A. Numerical model 



We consider the simplest model for our numerical simulations: d — 1 and we take a finite number of modes, with 
k^±n,l<n<N. We also take A(fc) = K{k) = 1 and R{k) = 1. We set = 1, = 1. For this choice of 
parameter our result (|43|) reads 



D, = D^- h 



^ ^' n—l 



I. — / 



We will perform four simulations: first we consider stochastic dynamics with detailed balance, and let the coupling 
h vary, to explore the range of validity of our perturbative result. Then, we will simulate the three special cases 
described above, for different values of the speed of evolution of the field n^. Each simulation is performed for 1 mode 
and 10 modes. 

In the simulations, we let one particle evolve in the field for a long time r ^ k^^, k^^, and we measure its position 
at a fixed set of times. We repeat this simulation a large number of times (around 10^) and, using these measurements, 
we compute {x{t)'^) /2dt, where t is the measurement time. For large t, this function fluctuates around a mean value, 
which gives us the effective diffusion constant. When these fluctuations are not small, they are taken into account 
with error-bars on the plots. 



B. Validity range of the perturbative result 



The first question that arises is: to what extent is our perturbative result valid? To find the validity range for h, we 
take stochastic dynamics with detailed balance, with 1 mode and — \ and look at the effective diffusion coefficient 
as a function of h. The comparison between the simulations and the result (j47p is given Fig. [1] and it shows that 
our computation is valid (i.e. the relative error is less than 5 %) for h < 1.2. A more relevant criterion is however to 
what extent the deviation from the bare result can be predicted by our result. Fig. [T] shows that the theory predicts 
the deviations of the diffusion constant from its bare value in the region where the diffusion constant deviates of the 
order of 15-20% from its bare value. 
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FIG. 1: Effective diffusion coefficient for stochastic dynamics witli detailed balance (Gaussian ferromagnet with model A 
dynamics) with a single mode field as a function of the coupling constant h. The crosses represent numerical simulations, the 
solid line is the perturbative result Eq. (|47[) . 




FIG. 2: Effective diffusion coefficient for stochastic dynamics with detailed balance (Gaussian ferromagnet with model A 
dynamics) as a function of with Dai = 1 and (a) 1 mode and h = 1 (b) 10 modes and h = 0.5: numerical simulations 
(crosses) and perturbative results (solid lines). 



C. Stochastic dynamics with detailed balance 

We can also vary the rate of evolution of the field k^. We set = 1 and plot Df^{K^) for 1 mode and h = 1, 
and for 10 modes and h = 0.5, comparing the numerical simulations results with (|47|) . The results are shown in Fig. 
[21 For 1 mode and /i = 1, at the border of the range of the perturbative approach, our results are in quite good 
agreement with the simulations. 
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(a) (b) 




FIG. 3: Effective diffusion coefficient for passive diffusion as a function of with Da, = 1 and (a) 1 mode and h — 1 (b) 10 
modes and h = 0.5: numerical simulations (crosses) and perturbative results (solid lines). 

D. Passive diffusion 

For passive diffusion, the results of numerical simulations are shown Fig. ^ and compared to Eq. ([Tf)) . We see that 
the analytical predictions are in good agreement with the results of simulations. In particular we see that depending 
on the relative values of and k^, the diffusion is either slowed down or speeded up. For small values of the 
diffusion is reduced, however on increasing the diffusion speeds up and passes through a maximum before decaying 
towards the bare value = 1 as predicted by our perturbative calculations. 

E. Particle not connected to any thermal bath 

The results for numerical simulations when the particle is not connected to any thermal bath are shown Fig. |4l 
They are in good agreement with Eq. (|47|) for > 1.5. When — >■ 0, according to our computation, the effective 
diffusion coefficient diverges, whereas physically it should go to zero, what is confirmed by the simulations. This 
discrepancy comes from the fact that we neglected the terms in o{tf) in our computation of {x{tf)'^Sint) and we took 
the limit — > oo before taking the limit — >■ 0. 

VI. DIFFUSION COEFFICIENT FROM THE PROBABILITY DENSITY FOR THE PASSIVE CASE 

Here we use a perturbation expansion of the Fokker-Planck equation in order to compute the perturbative correction 
for passive diffusion in a fluctuating field. The basic formalism is that the described in [H, 0] • In terms of the general 
model of this paper we are thus considering the special case C = and = = T . 

First, we give the elementary equations for the pure Brownian motion. The probability density function PQ{x,t) of 
a particle starting at a; = when t = 0, setting PQ{x,t < 0) = 0, satisfies 

Poix, t) = V • {D^VPix, t)) + Six)5{t). (48) 

We can write this equation in terms the free diffusion operator Hq — dt — V ■ (DxW): 

HoPo{x,t)=S{x)d{t). (49) 

and Fourier transform Po{k,uj) — J dx dtPi3[x^t)e^^^^'^^'^^^ is given by 



(50) 
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FIG. 4: Effective diffusion coefficient for a particle not connected to any thermal bath as a function of K0 with = 1 and (a) 
1 mode and h = 0.3 (b) 10 modes and h = 0.3: numerical simulations (crosses) and perturbative results (solid lines). 



The effective diffusion coefHcient of a process can thus be extracted from the Fourier transform of its probability 
density function P{k,uj) via P, Q 



-1 

1,2 1 



£) = (^^lim^fc^P(fc,0)j . (51) 

We are indeed interested in the large distance behaviour of our system, that is why the effective diffusion coefficient is 
given by the small wave-vector behavior. With this equation in hand, our strategy is simple: compute the probability 
density function for the passive diffusion and extract the effective diffusion coefficient. 

Now, we introduce a given field (j){x, t). For the general model described above, Eq. ([2]) gives the diffusion operator 
which replaces i/o hi Eq. P^ : 

H = Ha + H,,,t = dt -V ■ {D^V + hK^[V {K(l,){x,t)]). (52) 
The equation satisfied by the Fourier transform P{k,uj) is 

HP{k, u) = {iu; + D^e) P{k, io) + hn^ j ^^y+^ k ■ qK{q)4>{q, v)P{k -q,uj^u) = l, (53) 



so P(k,uj) is given by the integral equation: 



P{k,uj) = Po(fc,a;) 



J -^^^k- qk{q)4>{q,v)P{k- q,uj - v) 



(54) 



In this equation, the probability density function is that of pure Brownian motion, perturbed to the order h. Iterating 
this equation gives an explicit expression of P(k,uj) up to the desired order of h. 

Once we have the equation for P{k,uj) for a given field we need to extract the effective diffusion coefficient and 
proceed by averaging ([54]) over the configurations of the field (t){x,t) (which does not depend on the particle position). 
The field has a Gaussian probability density function with a two point correlation function that can easily be computed 
from Eq. ([5]): 

U{q, ^)4>{q', ^')) = . ^^^if.^j, X (2^Y+^6{q + q')5{v + ^')- (55) 
\ I + {K^R{q)A{q)y 

Now we have to insert this average in an explicit perturbative expansion of P{k, uj) given by (j54p . The lowest non- 
zero order is the second order: the field (f){x,t) has to appear at least twice. Moreover, we will obtain the propability 
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density function to the order h'^, what is exactly what we did with the path integral method. Computing P{k,uj) to 
the order and averaging the field out leads 



1 - 2h^Po{k,Lo)TK4,Kl 



dqdv k-q{k-q) qR{q)k{qf ~ 



{2nr+^ ^2 + («^i?^(q)A(q))2 



(56) 



Then, we can restrict ourselves to a; = 0, use the expression (|50p and integrate over v, using J 
1 



2tt (iiy+a)(i/2+/32) 



2/3(a+/3) ■ 



P{k,0)) = 



1 - 



dq 



k-q{k-q)-qk{qf 



{2nY A{q){DUk - q? + n^R{q)K{q)) 



(57) 



Finally, we just need to determine to behavior of the above expression when |fc| — >■ 0; a straightforward computation 
gives 



P(fc,0) 



1 



h^K^ r dq q^k{qf{D^q^ - K^R{q)K{q)) 



d J {2ttY K{q){D^q^ + K^R{q)K{qW 



(58) 



and we recover the effective diffusion coefficient given in (j45 



h Kif^ f d k 



r9'i^(9)% 



D^q^ - K^R{q)Aiq) 
A{q)[D^q^ + K^R{q)A{q)] 



(59) 



VII. CONCLUSIONS 



We have analyzed the diffusive behavior of a tracer particle diffusing in a time dependent Gaussian potential in the 
limit of weak coupling between the particle and the field. The method has the advantage that it can be applied to a 
wide range of models and not just the cases of passive diffusion or active diffusion with detailed balance. In these two 
aforementioned cases the method agrees with results obtained respectively via a perturbation expansion of the Fokker 
Planck equation and a Kubo formulation. We have also been able to look at non-equilibrium systems with variable 
feedback of the tracer on the field and also systems where the field and the tracer are subject to thermal noise of 
different temperatures. The range of behavior seen in the late time diffusion coefficient is quite rich, and depending 
on the models considered, coupling to the field can either slow down or speed up the diffusion. The speeding up or 
slowing down of diffusion and the possibility of a form of stochastic resonance depends on the relative rates of the 
dynamics of the fluctuating field and the bare diffusion constant of the tracer. Extensions of the work done here 
beyond the weak coupling approximation would be interesting to pursue, it is perhaps possible to apply Gaussian or 
mode coupling type approximations ?] to analyze the regime of strong interaction and perhaps even explore whether 
field fluctuations can lead to anomalous diffusion. In addition it would be interesting to see how the effects found 
here are modified when the coupling between the field and the tracer are non-linear, for instance quadratic. Such 
couplings are natural in systems where the tracer does not break the symmetry of the fluctuating field but rather 
enhances or suppresses its fluctuations. An example is a stiff membrane insertion which suppresses fluctuations in 
membrane curvature [ij] • A final point that would be interesting to address is what is the effect of a finite density of 
tracers for the active system under stochastic dynamics obeying detailed balance. As mentioned previously there will 
be induced interactions between the particles 23, 24] and it would be interesting to see how this modifies the effective 
diffusion constant. 
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